A Quantum-Theoretic Analog for a Pair of Noncommuting Observables of the 

Semiclassical Brillouin Function 



Paul B. Slater 

ISBER, University of California, Santa Barbara, CA 93106-2150 
e-mail: slater@itp.ucsb.edu, FAX: (805) 893-7995 
(February 1, 2008) 

We study, with the use of numerical integration, a noncommutative extension of a quantum- 
theoretic model (an alternative to the semiclassical Brillouin function) — recently presented by 
Brody and Hughston and, independently, Slater — for the thermodynamic behavior of a spin-i 
particle. Differences between the (broadly similar) predictions yielded by this extended model and 
those obtained from its conventional (semiclassical/ Jaynesian) entropy-maximization counterpart 
C ' are examined. 
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The Brillouin function, 

-£ , = tanh/3, (1) 
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ly-j ■ where E is the expected energy and [3, the inverse temperature parameter, has long served as a model of the ther- 
modynamic behavior of an ensemble of N noninteracting identical spin-i particles in an applied magnetic field [jp. 
In our simplified notation, we take (3 to represent where k is Boltzmann's constant, /i is the particle's magnetic 
moment, B is the external field strength, T is the temperature, and fiB — h, where h is Planck's constant. So, we set 

c^; h=1 - 

Recently, Brody and Hughston Q have argued that the theoretical underpinnings of ([[]) are semiclassical in nature, 
since the weighting of the phase space volume is eliminated, and random phases are averaged. Park and Band, in 
an extended series of papers questioned the conceptual foundations of the semiclassical/ Jaynesian approach to 
quantum statistical thermodynamics (cf. — from which (|l|) can be derived ||, p. 187]. Balian and Balasz || did, 
however, provide a rigorous justification (making use of field-theoretic concepts) for the Jaynesian (maximum-entropy) 
method, but, let us note that their argument was asymptotic in nature, relying upon a "supersystem" consisting of 
N replicas of the system, for which they required that N — > oo. 

Based on certain metrical considerations, Brody and Hughston proposed as a quantum-theoretic alternative to (Q), 
"pH . the function, 

~ E ~Wr (2) 

\ where the Ps represent modified (hyperbolic) Bessel functions. (This result was also presented — in a graphical 
manner — in a somewhat earlier paper of Slater Q.) Bessel functions often appear in the distribution of spherical 
and directional random variables ||. Ratios of modified Bessel functions, such as occur in (||), play "an important 
role in Bayesian analysis" [^J. It seems important to note, in this regard, the identity, 

tanh/3=-^ 7 -. (3) 

Lavenda j|, pp. 193 and 198] has argued, at considerable length, that the Brillouin function (0) lacks a suitable 
probabilistic basis, because the integral form for the modified Bessel function /„ exists only for v > ^. 

The relation (JsJ) has been used in expressing the expected energy of the linear-chain-lattice case (d = 1) of the 
Z?-vector (or n- vector) model for D — 1 [jl0[ p. 492], [O, p. 370], while the relation (||) emerges for the instance 
D = 4. The general expression in question takes the form, 

- E = 7^wy (4 » 

The reason that the specific dimension four, thus, arises in interpreting the analyses of Brody and Hughston O and 
Slater |Q, would appear to be due to the fact that the spin-i states (both mixed and pure) can be considered to lie 
on the surface of a hemisphere in /owr-dimensional Euclidean space, equipped with the (natural) Bures or statistical 
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distinguishability metric []i2|-|l4j . Let us also observe that the (classical) Langevin function 0] too is expressible as a 
ratio of modified Bessel functions, that is (.0 = 3), 

1 h(J3) 

coth e-~p = WY (5) 

We also point out that the analysis in j7| leads to the result (cf. (g)) corresponding to D = 6, 

for the /ive-dimensional convex set (the unit ball in five-space) of quaternionic two-level systems [ fl5f . (The result (Q) 
is based on the t/iree-dimensional convex set — the "Bloch sphere" , that is, the unit ball in three-space — of the 
standard complex two-level systems, which, as just mentioned above, also can be viewed as forming a hemisphere in 
four-space ]l2]|.) 

Brody and Hughston have argued that the differences (Fig. P in predictions between ([!]) and (||) might be tested 
in small systems (that is, small N), where "there seems to be no a priori reason for adopting the conventional mixed 
state approach." 




FIG. 1. The Brillouin function (|l|) and the (more steeply-sloped at /3 = 0) quantum-theoretic alternative 
between the two curves is of the greatest magnitude (.561292) at /3 — ±1.45489. 



The difference 



For an extended discussion of the role of negative temperatures, in this context, cf. fll| , sec. 3.52]. 

Brody and Hughston noted that the model (|2j) yielded a nonvanishing heat capacity at zero temperature. "Since it is 
known in the case of many bulk substances that the heat capacity vanishes as zero temperature is approached, it would 
be interesting to enquire if a single electron possesses a different behaviour, as indicated by our r esults" j|. They also 
observed that the increase in magnetization, when the temperature decreases, is slower for their quantum-theoretic 
result than for the scmiclassical one. 

In this letter, we extend the specific line of reasoning employed by Slater (?]] - based upon the Bures metric 
^ , ^6|Jl7| - to the case in which, rather than the expectation value (E) of one observable (as in (|J^1), one is 
interested in fitting the expectation values of two noncommuting observables (cf. |^Jl^|). We take these observables, 
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(7) 



to be two of the Pauli matrices. (One might also possibly use, Si = <Ti/2, S2 = 02/2, as the spin observables |lj], 
p. 38].) To obtain the conventional (semiclassical/Jaynesian) solution to this problem pp| , pl[ ], we express the target 
density matrix (p) as 



cxp(£! • / - Ai<ti - A 2 o- 2 ), 



(8) 



where Q + 1 and are the Lagrange multipliers for the normalization of p and the measured value of cti , respectively. 
These multipliers must satisfy 
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Q = — InTrcxp (— Ai<7i — A2CT2), 



(9) 



and 



{)X fa>- (* = 1,2) (10) 



The enforcement of these constraints leads to the result, 



Setting either Ai = or A2 = 0, we essentially recover the Brillouin function (jl]). 

Now, the volume element of the Bures metric over the three-dimensional convex set ("Bloch sphere") of spin-i 
systems is §§, eq. (6)] (cf. @||), 

(12) 



where 



— Hi -J- ^ 



is the additional Pauli matrix. If we integrate the term p2[ ) over one of the three coordinates (say, ((73)), we obtain 
simply a uniform distribution (7r/8) over the unit disk (0 < (01) 2 + (02) 2 < !)■ Interpreting this uniform distribution as 
a density-of-states or structure function, we can apply a bivariate Boltzmann factor, e -Pi( a i)-l 3 2( rT 2) ( to it. Integrating 

this product over the (<72)-coordinate (between the limits of ±\J 1 — (01} 2 ), we obtain, 

e-^^TTsinh (foyjl - (cti) 2 ) ^ 



4/3 2 

The corresponding partition function is the integral of ( |l4| ) over the remaining coordinate ((01}) from -1 to 1. This 
integration has to be performed numerically. Carrying this out, we are then able, again employing numerical inte- 
gration, to obtain the (twofold) expected value (((01}}) of (o~i) (Fig. ||) as a function of /3j and /?2, as well as the 
variance about this expected value (Fig. ||), and the covariance between (01) and ((72) (Fig. ||). (The covariance is 
the expected value — with respect to the Boltzmann distribution — of the product ((01) — ((01))) ((02) — ((02)))-) 
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FIG. 2. The expected value of the expected value of the observable <ti as a function of the inverse temperature parameters, 




FIG. 4. The difference between the quantum-theoretic results (Fig. y) and the semiclassical ones (Fig. |3|) for the expected 
value of the expected value of the observable ai — having identified the A's with the corresponding /3's 
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FIG. 5. The variance of the expected value of the observable <Ti as a function of the inverse temperature parameters, /3i and 
P2, of the quantum-theoretic model 




FIG. 6. The variance of the expected value of the observable o\ as a function of the inverse temperature parameters, Ai and 
A2, of the semiclassical (Brillouin-type) model 
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FIG. 8. The covariance between the expected values of the observables 01 and 02 as a function of the inverse temperature 
parameters, /3j and /?2, of the quantum-theoretic model 




FIG. 9. The covariance between the expected values of the observables o\ and <72 as a function of the inverse temperature 
parameters, Ai and A2, of the semiclassical (Brillouin-type) model 




FIG. 10. The difference between the quantum-theoretic results (Fig. |8j) and the semiclassical ones (Fig. |9|) for the covariance 
of the expected values of the observables o\ and o"2 — having identified the A's with the corresponding /3's 

For comparison purposes (cf. Fig. pi), we present the semiclassical (noncommuting Brillouin) counterparts to these 
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quantum-theoretic results in the companion figures (Figs. ||, g, . We note strong qualitative resemblances between 
the two members of each of these three pairs of figures. We also present in Figs. [|, [?] and [h], the differences obtained 
(after setting Ai to Pi, (i = 1,2)) by subtracting the semiclassical results (shown in Figs. |[ || and ^|) from the 
corresponding quantum-theoretic ones (given in Figs. ||, ^ and ||). The most substantial differences in all three cases 
appear in the vicinity of the (high-temperature) origin (/3i = (3% = 0) . 

Let us, in conclusion, consider the possibility of expanding the analyses above to the case of three noncommuting 
observables, rather than two. Then, we would apply a trivariate Boltzmann factor, e-^i^^-^i^^-P^^s) ^ ^ ^he 
volume element (|l^) itself of the Bures metric (rather than its two-dimensional [uniform] marginal — -|). Integrating 

out the (cr 3 )-coordinate (between the limits ±^/l — (ci) 2 — (02) 2 ), we obtain (cf. (jlj)), 

■ (15) 

where Jo is a Bessel function of the first kind. To obtain the corresponding partition function, it then appears 
necessary, similarly to before, to numerically integrate (|l^) over the unit disk (0 < (ax) + (a 2 ) < 1). 
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